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We report an analysis ol charmless hadronic decays of neutral B mesons to the final state 
KsK ± ir T , using a data sample of (465 ± 5) x 10 6 BB events collected with the BABAR detec- 
tor at the T(4S) resonance. We observe an excess of signal events with a significance of 5.2 
standard deviations including systematic uncertainties and measure the branching fraction to be 
B (B° — > KsK ± ir^ 1 ) = (3.2 ± 0.5 ± 0.3) x 10 -6 , where the uncertainties are statistical and system- 
atic, respectively. 

PACS numbers: 13.25.Hw, 14.40.Nd 



Charmless decays of B mesons to hadronic final states tions. In recent years, there has been a surge of new 
containing an even number of kaons are suppressed in results on these decays: B° — > K®K® and B + — > K®K + 
the standard model (SM). Decays of this type mainly have been observed [H, E|, and there is evidence for the 
proceed via the b — > d "penguin" transition, involving related vector-vector final states Only upper limits 
a virtual loop, and hence are sensitive to potential new on the corresponding pseudoscalar-vector final states ex- 
physics contributions since the presence of new particles ist: B(B° K °K *°) +B(B° -> K°K*°) < 1.9 x 10~ 6 [1] 
in the loops can produce deviations from SM expecta- and B(B + K + K*°) < 1.1 x 10~ 6 [3], both at 90 % con- 
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fidence level (unless explicitly stated otherwise we use the 
symbol K* to denote the if* (892) resonance and the in- 
clusion of charge conjugate modes is implied). Note that 
decays with additional suppression in the SM, such as 
B° -> #(*)+#(*)-, which are expected to proceed via 
annihilation amplitudes, have not been observed 

m. 

Since the vector resonances involved have non- 
negligible widths, the pseudoscalar- vector decays are best 
studied using Dalitz plots of the three-body KKir final 
states. In the three-body channels, contributions from 
suppressed b — > u tree amplitudes are expected to be 
important, in addition to the b — > d penguin amplitudes. 
Recent investigations of three-body channels suggest that 
additional resonances are present. Most notably, the 
B + — > K + K~n + channel exhibits an unexpected peak 
near 1.5GeV/c 2 in the K + K~ invariant- mass spectrum, 
which accounts for approximately half of the total event 
rate We call this peak, with unknown spin and 

isospin quantum numbers, the /x(1500). The lack of a 
/x(1500) signal in B + — > K®K®ir + decays implies that 
the /x(1500) does not have even spin if isospin is con- 
served in the decay (ljj . A search for an isospin partner 
to the /x(1500) that decays to K°K + and which could 
be produced recoiling against a pion in B decay could 
help to clarify the nature of this resonance. 

In this paper, we present the results of a search for 
the three-body decay B° — > K^K^ir^, including inter- 
mediate two-body modes that decay to this final state 
but do not contain charm quarks. No decays to this 
final state have been observed as yet. The best avail- 
able upper limit on the inclusive branching fraction is 
B(B° ->■ K^K^if?) < 18 x 1(T 6 0. There appears to 
be no explicit prediction for the inclusive branching frac- 
tion of J5° — > K^K^tt^ . Some theoretical predictions ex- 
ist, however, for the relevant resonant modes. Expected 
branching fractions for B° -> (K*°K° + K*°K a ) and 
B° -> K*±K* are in the range (0.2-2.0) x 10~ 6 and 
(0.2-1.0) x 10 -7 , respectively |17H24| . Extensions to the 
SM can yield significantly larger branching fractions. For 
instance, in supersymmetric models with i?-parity viola- 
tion, the branching fraction for B° -> (K*°K° + K*°K°) 
could be as large as 10~ 5 [Hj]. 

The data used in the analysis, collected with the 
BABAR detector [26[ at the PEP-II asymmetric energy 
e + e~ collider at SLAC, consist of an integrated luminos- 
ity of 424fb _1 recorded at the T(4S) resonance ("on- 
peak") and 44fb _1 collected 40 MeV below the reso- 
nance ("off-peak"). The on-peak data sample contains 
(465 ± 5) x 10 6 BB events. 

We reconstruct B° — > K^K^n^ decay candidates by 
combining a K® candidate with one charged kaon and 
one oppositely charged pion candidate. The and 
tt^ candidates are required to have a minimum trans- 
verse momentum of 50 MeV/ c and to be consistent with 
having originated from the interaction region. Identifi- 



cation of charged kaons and pions is accomplished with 
energy-loss information from the tracking subdetectors, 
and the Cherenkov angle and number of photons mea- 
sured by a ring-imaging Cherenkov detector. We dis- 
tinguish kaons from pions by applying criteria to the 
product of the likelihood ratios determined from these 
individual measurements. The efficiency for kaon selec- 
tion is approximately 80 % including geometrical accep- 
tance, while the probability of misidcntification of pions 
as kaons is below 5 % up to a laboratory momentum of 
4 GeV/c. A K® — > tt + tt~ candidate is formed from a pair 
of oppositely charged tracks (with the pion mass hypoth- 
esis assumed) having an invariant mass that lies within 
15 MeV/ c 2 of the nominal Kg mass 27], corresponding 
to 5 times the Kg mass resolution. We require the ratio 
of the measured Kg decay length and its uncertainty to 
be greater than 20, the cosine of the angle between the 
line connecting the B and K® decay vertices and the Kg 
momentum vector to be greater than 0.999, and the Kg 
vertex fit probability to be greater than 10~ 6 . 

To suppress the dominant background contribution, 
which arises from continuum e + e~ — > qq (q = u, d, s, c) 
events, we employ a Fisher discriminant that combines 
four variables. These are the ratio of the second to the 



zeroth order momentum- weighted angular moment [28j, 
the absolute value of the cosine of the angle between the 
B direction and the beam axis, the magnitude of the 
cosine of the angle between the B thrust axis and the 
beam axis, and the proper time difference between the 
decays of the two B mesons divided by its statistical un- 
certainty. The first three quantities are calculated in the 
center-of-mass (CM) frame. 

In addition to the Fisher output (J-), we distinguish 
signal from background events using two kinematic vari- 
ables: the difference AE between the CM energy of the 
B candidate and y/s/2, and the beam-energy-substituted 
mass ities = \/ s /4 — P^, where y/s is the total CM en- 
ergy and pb is the momentum of the candidate B meson 
in the CM frame. The signal tties distribution peaks near 
the B mass with a resolution of about 2.6MeV/c 2 , while 
its AE distribution peaks at zero with a resolution of 
approximately 20 MeV. We select signal candidates that 
satisfy 5.272 < m ES < 5.286 GeV/c 2 , \AE\ < 0.075 GeV, 
and T > —0.145. The requirement on T removes approx- 
imately 70 % of continuum background while retaining 
90 % of signal events. 

Another source of background arises from B decays, 
mostly involving intermediate charm or charmonium 
mesons, or charmless final states that are misrecon- 
structed. We exclude B candidates that have two-body 
mass combinations in any of the following invariant-mass 
ranges: 1.82 < miK^R*) < 2.04, 1.81 < m(K°n^) < 



1.91, 1.83 < m{K ± TT^) < 1.90, 3.06 < to^tt^) < 
3.17, and 3.66 < m(K ± 7T^) < 3.73 (all in units of 
GeV/c 2 ). These ranges reject decays from D + and Df, 
D + . D°, J/tp, and i/>(2*S f ) mesons, respectively. Charmo- 
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nium contributions result mainly from the leptonic de- 
cays of J/tfj and i/)(2S), where one lepton is misidentified 
as a charged pion and the other as a kaon. 

The efficiency for signal events to pass all the selec- 
tion criteria is determined as a function of position in 
the Dalitz plot. Using a Monte Carlo (MC) simulation 
in which events uniformly populate the phase-space, we 
obtain an average efficiency of 20%, though values as 
high as double (as low as half) that value are found near 
the center (corners) of the Dalitz plot. 

An average of 1.1 B candidates is found per selected 
event. In events with multiple candidates we choose the 
one with the highest B vertex fit probability. We verify 
that this procedure does not bias our fit variables. In 
some signal events, the B candidate is misreconstructed 
due to one track being replaced with a track from the 
rest of the event. The fraction of such events is below 
2% in the phase-space MC, but is closer to 5% in MC 
samples where the events populate the K* bands. Mis- 
reconstructed signal events are considered as a part of 
the signal component in the fit described below. We as- 
sign a systematic error to account for the uncertainty in 
the rate of these events, which is related to the unknown 
Dalitz-plot distribution of the B° — > K®K ± n T decay. 

We study residual background contributions from BB 
events that survive the invariant-mass exclusion require- 
ments described earlier, using MC simulations. It is 
found that these events can be combined into four cat- 
egories based on their shapes in m,Es and AE. The 
first category {BBi) comprises B° — > rj'K®,r]' 
and misreconstructed B° —> D~tt + ,D~ 
cays and has a broad peak in toes and a nonpeaking 
AE shape. The second and third categories (BB 2 and 
BB$) represent the charmless decays B a —} K®K + K~ 
and B° — > K®ir + ir~ , where a kaon or a pion is misidenti- 
fied leading to a AE distribution that peaks with nega- 
tive or positive mean, respectively. The MC simulations 
of these decays are based on our recent studies of their 
Dalitz plot distributions [29|, |3(j. The fourth category 
(BB^) contains the remainder of the BB background 
and is mainly combinatorial in nature. Based on the 
MC-derived efficiencies, total number of BB events, and 
known branching fractions 27, |3l|, we expect 25, 173, 
215, and 668 events from the four BB background cate- 
gories, respectively. 

To obtain the B° — > K^K ± tt t signal yield, we per- 
form an unbinned extended maximum likelihood fit to 
the candidate events using three input variables: toes, 
AE, and T . These variables are found to be largely 
uncorrelated - the maximum correlation is between the 
signal toes and AE distributions and is about 13 %. For 
each component j (signal, qq background, and the four 
BB background categories), we define a probability den- 
sity function (PDF) 



where i denotes the event index. The extended likelihood 
function is given as 



n 

K° S K~ de- 



E 



n-P) 



(2) 



where rij(fc) is the yield of the event category j{k). 

For the signal component, the toes and AE distribu- 
tions are each parametrized by the sum of two Gaussian 
functions, while the J- distribution is described by a bi- 
furcated Gaussian function with a small admixture from 
the sum of two Gaussians. We fix the shape parame- 
ters to the values obtained from the B° — > K^K ± tt t 
phase-space MC sample, after adjusting them to account 
for possible differences between data and MC simulations 
determined with a control sample of B Q — > D~tt + , D~ — > 
A"!?7r~ decays. For the continuum background, we use an 



ARGUS function 32j to parametrize the toes shape and 
a linear function for AE. The continuum Fisher shape 
is modeled with a function that is composed of a Gaus- 
sian tail with relative fraction 99.6% (large component) 
and a small Gaussian with different mean and width val- 
ues. This shape provides a good description of the off- 
peak Fisher distribution, as well as of the corresponding 
MC distribution. One-dimensional histograms are used 
as nonparametric PDFs to represent all three fit variables 
for the four BB background components. 

The free parameters of our fit are the yields of signal, 
BB2, BB%, and continuum background together with 
the slope of the continuum AE PDF and the mean and 
width of the large Gaussian component of the continuum 
T PDF. The ARGUS £ parameter and parameters of 
the small Gaussian component of the continuum Fisher 
function are fixed to values determined from candidates 
selected in the off-peak data sample with a looser re- 
quirement on toes- The yields of BB\ and BB4, and all 
shape parameters of the four BB background categories 
are fixed to the values determined from MC simulations. 

We cross check our analysis procedure by removing 
the requirements that reject backgrounds from B de- 
cays involving charm mesons, instead selecting regions 
of the Dalitz plot dominated by intermediate charm 



mesons. We select B' 



D~ 



,D- 



K°K~ de- 



V i j = V j {mvs i )'PA&E i )V j {r 



(1) 



cays requiring 1.84 < miK^K^ < 1.89GeV/c 2 , and 
B° D~K + ,D~ K°ir~ decays requiring 1.84 < 
TO( J ftr°7r ± ) < 1.89GeV/c 2 . We then apply our fit to find 
the yields for the B° ->■ D--K+ and B° ->■ D~ K+ chan- 
nels. We find values consistent with the expectations 
based on world-average product branching fractions [27I ] 
within statistical uncertainties. 

Applying the fit method described above to the 14 276 
selected candidate B° — » K°K ± tt t events, we find 
262 ± 47 signal events. The fitted yields of the BB 2 
and BB3 categories are 199 ± 51 and 262 ± 55, respec- 
tively, consistent with the MC-based expectations. The 
fitted values of all other free parameters of the fit are also 
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consistent with expectations based on studies of control 
samples and MC simulations. The results of the fit are 
shown in Fig. [TJ The statistical significance of the signal 
yield, given by the square root of the difference between 
twice the value of negative log likelihood obtained assum- 
ing zero signal events to that at its minimum, is 6.0 a. 
Including systematic uncertainties (discussed below) , the 
significance is 5.2 a. 

The B° — > K^K^tt 11 branching fraction is determined 
from the result of the fit by calculating signal proba- 
bilities for each candidate event with the s Vlot tech- 
nique [33(. These are divided by event-by-event efficien- 
cies that take the Dalitz-plot position dependence into 
account, and summed to obtain an efficiency-corrected 
signal yield of 1326 ± 207 events. We further correct 
for the effect of the charm and charmonium vetoes (es- 
timated using a range of MC samples with different 
Dalitz-plot distributions), and divide by the total num- 
ber of BB events in the data sample assuming equal 
production of B°B° and B + B~ at the r(4S). The re- 
sult for the branching fraction is B (B° — > K®K ± 7r =F ) = 
(3.2 ± 0.5 ± 0.3) x 10~ 6 , where the first error is statistical 
and the second is systematic. 

We find the systematic error to be due to uncertainties 
in the signal PDFs (5.2%), including possible data-MC 
differences in the signal PDF shapes evaluated using the 
control sample of B° —> D~n + ,D~ —> K®tt~ decays; 
uncertainties in the background PDFs (2.5 %), including 
effects due to the fixed values of some of the qq PDF 
parameters (recall that the parametrization used is vali- 
dated with off-peak and MC samples and that the most 
critical parameters are floated in the fit to data; the un- 
certainties are evaluated by varying the fixed parame- 
ters) and due to the fixed content of the histograms used 
to describe the BB background PDFs; potential fit bi- 
ases, studied using ensembles of simulated experiments 
where continuum events are drawn from the PDF shapes 
and signal and BB background events are randomly ex- 
tracted from MC samples (1.1%); uncertainties in the 
efficiency due to tracking (0.8%), K° s selection (0.9%), 
and particle identification (2.8%); and the error in the 
number of BB events (1.1%). We assign two system- 
atic uncertainties to account for the nonuniform Dalitz 
plot structure of the signal, both of which are estimated 
from MC simulations with different resonant contribu- 
tions: uncertainty in the fraction of misreconstructed 
events (3.0%) and uncertainty in the correction due to 
vetoes (4.1%). Other sources of systematic uncertainty, 
including the fixed yields of BB\ and BB4, are found 
to be negligible (recall that the fitted yields of BB2 and 
BB 3 are consistent with expectation). 

In Fig. [2] we show the efficiency-corrected Dalitz plot 
for signal decays, obtained using event-by-event signal 
probabilities. We verify that this technique correctly re- 
constructs the signal Dalitz plot distribution using MC 
simulations in which the B° — > KgK^Tt^ events contain 



different structures. There appears to be some structure 
in the K*° region at low J ftT ± 7r =F invariant mass, and an 
excess of events at low K^K^ mass with a highly asym- 
metric hclicity angle distribution. Quantitative state- 
ments concerning the content of the Dalitz plot require a 
dedicated amplitude analysis, which is beyond the scope 
of the present study. However it appears that there is 
no major contribution from an isospin partner of the 
/x(1500) decaying to K®K + , which contrasts to the clear 
signal seen in B + — » K + K~n + decays 14 1. 

In summary, using the full BABAR data sample of 
424 ftT 1 collected at the T(4S) resonance, we find evi- 
dence for charmless hadronic decays of neutral B mesons 
to the final state K^K^n^. The signal has a sig- 
nificance of 5.2cr, after taking systematic effects into 
account. We measure the branching fraction to be 
B (B° -> K^K^ir?) = (3.2 ± 0.5 ± 0.3) x 10" 6 . We con- 
vert this result to B[B° (K°K+Tr- + K°K--k + )] = 
(6.4 ± 1.0 ± 0.6) x 10~ 6 by multiplying by a factor of 2. 
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FIG. 2: Efficiency-corrected Dalitz plot distribution of B — > 
K^K^ix^ decays, obtained with the s Vlot technique [33l ]. 
Bins with negative content appear empty as do regions corre- 
sponding to the charm and charmonium vetoes. The area of 
each box is proportional to the number of weighted events in 
that bin; the largest bin corresponds to 12 events. 
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